A tunable coupler with ScS quantum point contact 
to mediate strong interaction between flux qubits 



o 

(N 
(N 



A. A. Soroka 1 and V.I. Shnyrkov 2 

National Science Center "Kharkov Institute of Physics and Technology", 
Akhiezer Institute for Theoretical Physics, 61108 Kharkov, Ukraine 
2 B. Verkin Institute for Low Temperature Physics and Engineering, 
National Academy of Sciences of Ukraine, 61103 Kharkov, Ukraine 

In this paper we propose a kind of quantum inductance couplers (QUINC) which represents 
a superconducting loop closed by ScS quantum point contact, operating in deep quantum low- 
temperature regime to provide tunable (Ising-type) ZZ interaction between flux qubits. This coupler 
is shown to be well tunable by an external control magnetic flux and to provide large inter-qubit 
interaction energies | J/fce | — 1 K thus being very promising as a qubit-coupling device in a quantum 
register as well as for studying fundamental low-temperature quantum phenomena. Some entangle- 
ment measures of a two-qubit system are analyzed as functions of inter-qubit interaction strength. 

PACS numbers: 03.67.Bg, 74.50.+r, 85.25.Cp 
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I. INTRODUCTION 



Handling interaction between basic elements (qubits) 
of a quantum computer is one of key problems for im- 
plementation of computation algorithms thereinir— . In 
a computer register including N qubits, computation 
operations are generally unitary transforms U(2 N x 2 N ) 
of a superpositional state vector in the 2 Ar -dimensional 
tensor-product Hilbert state space of all qubits with the 
basis: {si} ® {s 2 } ® • • • ® {sjv}, {st} = {|t),- ! \i) J be- 
ing basis sets of separate qubits. Unitary transforms 
U (2™ x 2 n ) of a quantum register state vector in 2™- 
dimensional subspaces which are produced by turning on 
interaction within certain subsystem of n qubits during 
some timespan and which realize specific computation 
operations are called n-qubit gates. A universal set of 
gates for quantum computation is such that generates all 
possible unitary transformations in the full Hilbert vector 
space and thus suffices for implementation of an arbitrary 
algorithm. In quantum informatics, the Brylinski's theo- 
rem states that a universal set of computation gates may 
be constituted from all one-qubit gates (providing local 
unitary transformations of separate qubits) and any non- 
primitive, that is entangling, two-qubit gate^. The most 
known entangling two-qubit gate in quantum informatics 
is the CNOT gate 5 . Different entangling two-qubit gates 
are interrelated involving one-qubit gates^. 

So, a physical system used for building a quantum reg- 
ister should provide a tuning of interaction energy be- 
tween pairs of qubits to form regulated entangled two- 
qubit states, with the possibility of turning on/off inter- 
qubit interaction^. A coupler with a quantum point con- 
tact (QPC) we propose realizes ZZ two-qubit tunable in- 
teraction generating a class of entangling gates that, as 
follows from the aforesaid, produce a universal set for 
quantum computation. Furthermore, control of interac- 
tion between the quantum coherent systems is of great 
interest to study fundamental effects of entanglement and 
correlations in the states of coupled qubits and qutrits 1 . 



Lately, a system of two coupled Josephson, and par- 
ticularly, flux qubits has been intensively studied 8 . In 
first papers simple systems of flux qubits were studied, 
with qubits located side-by-side and coupled by constant 
coefficients of mutual inductances My, determined by 
geometry of the relative position and self-inductances of 
qubits 9,10 . At the same time, new versions of coupling 
elements were proposed theoreticall y 11 ' 12 which, based 
on dc- and rf- SQUIDs, enabled to vary the magnitude 
and the sign of magnetic interaction of flux qubits to 
be coupled by means of an external control parameter 
(the bias current in dc SQUID and the external mag- 
netic flux applied to the loop of rf SQUID). These sign- 
and magnitude-tunable qubit-coupling elements were im- 
plemented in the Refs. [l3lflil for the first time, but ob- 
tained magnitudes of qubits' interaction energies, pro- 
portional to respective dynamical Josephson inductances, 
were rather small (about 50 mK), being insufficient for 
coupling qubits with high tunnel energy splitting^. 

In this paper we propose a tunable coupler to provide 
strong ZZ interaction between flux qubits using the quan- 
tum inductance of a superconducting loop closed by ScS 
quantum point contact^. A quantum inductance coupler 
(QUINC) is a nonlinear quantum system being in some of 
its eigenstates, such that additional unwanted entangle- 
ment between this state and states of the coupled qubits 
could be excluded. On this account, the coupler working 
state, with largely varying quantum inductance (curva- 
ture of the energy vs. flux dependence), is found that 
forms in a three-well asymmetric potential of the quan- 
tum loop in such a way as to be localized in a side well of 
the potential and tunnel to its central well very weakly. 
We study the quantum inductance of such coupler, which 
determines interaction energy between two flux qubits, 
depending on the external magnetic flux through the cou- 
pler's loop as the control parameter. Also we analyze 
some widely-used entanglement measures for a system of 
two identical qubits with tunable coupling energy of ZZ 
type as functions of its magnitude. 
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II. MODEL AND NUMERICAL ANALYSIS 



Let's consider a system of two flux qubits having tun- 
nel splittings AE^' = 2Ai and AE^ — 2A 2 in points of 
symmetry of their two-well potentials, which are induc- 
tively coupled to the coupler via transformers of magnetic 
flux with coefficients of mutual inductance Mi and M 2 , 
respectively. The coupler we select is a superconducting 
Nb 3D-loop closed by a clean ScS quantum (atomic-size) 
point contact, such that d~ A^ , d<ti I, £o, where d is the 
contact dimension, is the electron wave length, I is the 
electron elastic mean free path, £ is the superconducting 
coherence lengt h 17 ' 18 . Such 3D-loop of the coupler may 
be topologically the same as that we used for designing 
the qutrhJ£, but with a different set of parameters: the 
loop inductance L, the self-capacitance C and the criti- 
cal current I c of the ScS contact (L and I c determining 
a value of the nonlinearity parameter /3l = 2irLI c /$o, 
$o = h/2e is the flux quantum), that is necessary to es- 
tablish peculiar properties of the coupler's energy struc- 
ture. Under operation conditions, the coupler is initial- 
ized by means of external flux <f> e to certain energy level 
E n (<& e ), which nonlinear dependence on control parame- 
ter $ e determines coup ling strength between qubits. 

It was shown in Ref.[l2|, that energy of Ising-type ZZ 
interaction of two flux qubits, mediated by a coupler as 
a nonlinear quantum element, is 



J($ e ) = -x($e)M 1 M 2 I 1 I 2 , 

J_ dl n _ d 2 E n _ 1 
Lq d* fi d$? L 



1 



(1) 



where I\ , 1% are the superconducting currents of the ba- 
sis states {|t) j j I4-)*} (i = 1 , 2) of the first and the second 
qubits; _E„($ e ) and /„($ e ) are the energy of the cou- 
pler operation level and circulating in its loop supercon- 
ducting current as functions of magnetic flux $ e applied 
to the loop, respectively; $($ e ) = <& e — LI n (<&) is the 
internal magnetic flux in the loop; x($ e ) = E'^($ e ) is 
the local curvature of the coupler operation level called 
the susceptibility or, equivalently, the reciprocal quan- 
tum inductance Lq 1 of the level. As seen from Eq. ([1}, 
interaction energy </($ e ) between two flux qubits is pro- 
portional to the curvature function x(<I> e ). Positive val- 
ues of the curvature (x > 0) correspond to ferromagnetic 
coupling of qubits (FM, J < 0), while negative values 
(x < 0) to antiferromagnetic coupling (AFM, J > 0). It 
should be emphasized that a quantum inductance cou- 
pler, being a quantum element, must function in one- 
dimensional Hilbert space, so as not to get entangled with 
the qubits, thus resulting in effective four-dimensional 
Hilbert state space of the qubit subsystem. It can be 
achieved provided that (i) the operation level E n ($ e ) 
is weakly-superpositional, (ii) the distances between E n 
and the neighboring levels E n ^i, E n+ i substantially ex- 
ceed the splittings AEq 1 ^ , AEq± of the qubits. Note that 
quantum coherent adiabatic regime of the coupler func- 
tioning allows to eliminate decohering influence of quasi- 
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Figure 1: The coupler potential U (f)/ks (in units of tempera- 
ture) at f3 L = 4.4, C = 3.1 fF, f e = 0.018; E„ is the operation 
level of the coupler, | V&n. | 2 is the squared wave function of the 
respective quantum state. 



particle currents, inherent to a classical system (SQUID), 
on the qubit dynamics. 

To analyze properties of the QUINC with ScS atomic- 
size contact (such that the parameter g = Ej/Ec — 
&oIcC / (2ire 2 ) ^> 1) at the bath temperature T con- 
siderably less than any distances between adjacent en- 
ergy levels of the quantum system (thus validating 
a zero-temperature approximation), we use the flux- 
representation Hamiltonian in the form 151 " 19 



P 2 



Hc=^ + U(f;f e ) = 



h 2 d 2 

2M~dp 



$qIc 
2tt 



-2|cos(tt/)| 



27r 2 (/-/ e 



PL 

(2) 

where / = ( E ) /<I ) o and f e — ®e/&o are the normalized 
internal magnetic flux $ in the coupler loop and exter- 
nal magnetic flux $ e applied to the loop, M = ^ 2 C is 
the respective effective mass. The quantum dynamical 
observable of the internal magnetic flux in the coupler 
loop is given by an operator $ conjugated to an operator 
Q of charge in the contact capacitance: [$, Q] = ih^. 
The key feature of Hamiltonian ([2]) is its singular poten- 
tial U{f; f e ) following from the non-sine current-phase 
relation of ScS quantum point contact^ 8 -: 

I a (tp) = J c sin(<p/2)sgn[cos(<p/2)], I c = = N^-, 

eitN n 

where Ao is the superconducting energy gap, Rn is 
the normal-state resistance of the contact, N is an in- 
teger. The critical current of ScS atomic-size contact 
is quantized in consequence of quantization of the con- 
tact conductance i?^ 1 in units of Go = 2e 2 /h, that was 
observed experimentally 17 . Note that such an effective 
quantum Hamiltonian @ of the superconducting loop 
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Figure 2: The function (LL Q 1 )(f e ) for L = 0.3nH, fa =4.4, 
and the values of C: 1 - 3.0 fF, 2 - 3.05 fF, 3 - 3.10 fF. 

closed by the ScS contact, with its Josephson energy 
(Uj((f) = — (I c $o/tt)\ cos(v?/2)|) having the singular pe- 
culiarity and the dissipation vanishing at zero tempera- 
ture, satisfactorily describes the experiments both (i) on 
macroscopic quantum tunneling phenomena in the loop 
with the clean ScS contac t 20 ' 21 and (ii) on coherent quan- 
tum superposition of macroscopically distinct states in 
the three-well potential of the superconducting qutriti£. 
The solutions of the stationary Schrodinger equation 

H C (f;fe)W) = E(f e )*(f) (3) 

with Hamiltonian @ describe wave functions and 
energies E(f e ) of the stationary states of the coupler at 
a specified external magnetic flux f e . Let us consider 
solutions of Eq. ([3]) with the following coupler parame- 
ters: L = 0.3nH, C = 3.1fF, j3 L = 4.4 (I c «4.8//A) and 
f e = 0.018, corresponding to asymmetric three-well po- 
tential U(f;f e ), shown in Fig.[TJ In this configuration, 
we are interested in the state with minimal energy level 
E n lying both in the right-side and the central wells of 
the potential. Experimentally, this state ^ n (f) may be 
obtained from the initial ground state in the three-well 
symmetric potential at 3> e = $ as a result of decreasing 
$ e to the required value. With decreasing $ e , a central 
potential well transforms into a side well, where the qua- 
siclassical state with a large quantum number (n = 25 for 
E n in Fig.[IJ is localized. The level E n lies far below the 
potential barrier top, (Ub — E n )/kB ~ 12.7 K, and sep- 
arated from neighboring energy levels by a considerable 
interval AE: (E n - £ B -i)/*B « (E n +i - E n )/k B a 6K. 
As the state ^f n (f) is almost entirely localized in the 
right well, it proves to be utterly weakly-superpositional. 
For the parameters in Fig.Q] [such that x(fe) ~ 0] the 
probability of the state ^ n (f) being in the right well 
Pr = J R \^n(f)\ 2 df « 0.97, and in the central well 
P c = f c \t> n (f)\ 2 df w 0.03. The quantity P c charac- 
terizes a superpositional leakage of the state 4 r n (/) to 
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Figure 3: (a) The function (LLq 1 )^) for L = 0.3nH, C = 
3.10 fF and the values of fa: 1 - 4.4, 2 - 4.2; (b) The function 
-J(f e )/k B for L = 0.3nH, C = 3.10fF, Mi = M 2 = 0.03 nH, 
h = h = 0.8fj,A, and the values of fa: 1 - 4.4, 2 - 4.2. 



the central well. Further the coupler is characterized in 
an operation / e -range which is defined by a criterion of 
ensuring a desired fixed range of its normed susceptibil- 
ity [(LLg X X/ e ) € (-10 . . 10), see below]. When varying 
f e away from the point fo of zeroing of the ^-function 
(at a given /3l), the Pc magnitude increases; for curves 
in Fig. [2 the maximal Pc — 0.04 in the operation f e - 
range. With increase in 0l the potential barrier between 
the right and the central wells grows (and so grows the 
depth of the right well), that causes reducing of the Pc; 
e.g. for (3l — 4.6 the maximal Pc = 0.02 in the operation 
range. At that the Pr approaches nearer to unity. 

It should be pointed out that under actual conditions 
the coupler state will be always to some extent deco- 
hered (being in the so-called partially-coherent state), 
because of finite nonzero environment temperature and 
influence of unavoidable decohereing (noise-originating) 
factors. We assume the situation in which the influence of 
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small superpositional probability Pc of the coupler being 
in the central well on its entanglement with the qubits 
is smoothed away due to the finite noise-induced phase 
dispersion of the wave function ^„(/). Thus, weak super- 
positionality of the coupler operation state and its large 
separation from the neighboring states (multiply exceed- 
ing the tunnel splittings of qubits) can make the coupler 
a well-defined entity not entangled with the qubits. 

Because of small setting times of the operation state 
during variation of the external magnetic flux, estimated 
as r~ h/AE ~ 10 -11 s, the coupler will be described by 
the adiabatic susceptibility x($ e ) even at rather high 
switching rates. Similarly to stability of the base super- 
position state of the qutrit toward relaxation to underly- 
ing states 1 ^, the coupler operation level will be the same 
stable under certain conditions. It is possible due to de- 
signing of the coupler loop in the form of a high-quality 
three-dimensional toroidal superconducting cavity with 
no resonant modes at frequencies corresponding to that 
of transitions from the operation level to underlying ones. 

For practical analysis, the function (LLq )(/ e ) of the 
susceptibility (i.e. Lq 1 ) normed to the loop inductance 
L is convenient 16 . It is obtained by numerical solving 
Eq. ([3]), finding E n (f e ) and its double differentiation by 
the f e parameter: (LLQ\f e ) = (L/<S> 2 )d 2 E n /df 2 . InFig.H 

the function (LLq )(/ e ) is shown for /3l — 4.4 and sev- 
eral values of the ScS-contact capacitance. With increas- 
ing C, the curve [LLq )(/ e ) shifts right parallel to the 
x-axis. At that a form of the curve does not change 
qualitatively in the operation / e -range of the coupler. In 
Fig-ISJa) the function (LL^X/e) is shown for two differ- 
ent values of (3l- As discussed above, with increasing 
Pl the state *ff n (f) still more localizes in the right well. 
Nonlinearity of the E n (f e ) function increases, and so the 
(LLg^/e) rises in absolute value, with its slope at the x- 
axis decreased. This decrease of the susceptibility slope 
about the /o point can be compensated by widening of 
the coupler operation range as necessary. Fig. [2Kb) shows 
the function — J(/ e )/fcB (see Eq. Q])) of the mediated by 
the coupler inter-qubit interaction energy (in tempera- 
ture units and with minus sign for convenience) for the 
same set of the coupler parameters as in Fig.[3Ia) and the 
characteristic absolute values of coupler-qubit mutual in- 
ductances Mi = M 2 = 0.1L and of the basis currents in 
qubits. It is seen that at chosen parameters the variation 
of the control external flux in the relatively small range 
f e G (0.014 . . 0.022) ensures the wide range of variation 
(— J//cb) € (— 1 . . 1 K) of the inter-qubit coupling energy. 



III. ENTANGLEMENT MEASURES 

Let us now analyze properties of the two-qubit sys- 
tem with ZZ interaction J(/ e ) which may be provided 
by the considered coupler. Having obtained eigenstates 
and density matrix of the system, we will find explicitly 
entanglement measures widely used for characterization 



of the entangled states^ 2 -—. The effective Hamiltonian 
of the system in the tensor-product basis 

{It)! ; |4) 1 }®{|t> 2 ; ll> 2 } = (Iff) ; in) ; lit) ; in)} (4) 

(the numbers of arrows correspond to omitted qubit in- 
dices in " two-arrow" states) of the four-dimensional state 
space, close to qubits' degeneracy points, have the form 

#12 = (-£i<5>-Aio- x ) ® (7 + (T ® (-e 2 a z -A 2 a x ) + 

+ J(fe) O z ®G z = 



-62- 


-J -A 2 


-Ai 





-A 2 


-e\+S2- J 





-Ar 


-Ai 





£l-£2- J 


-A 2 





-Ai 


-A 2 £l - 


+-£2- 



(5) 

Here E\ , £2 are the energy biases of first and second qubits 
relative to their degenerate energy levels in symmetric 
two- well qubit potentials, <j Xi <j Zi oq are the Pauli matri- 
ces. The fourth-order characteristic equation for finding 
eigenvalues of H\ 2 [Eq.©] 

E 4 - 2(e 2 1 +e 2 2 + A 2 1 + A 2 2 + J 2 )E 2 - (8 £l e 2 J)E 
+ [e\+e 4 2 + J 4 + (A 2 l -A 2 2 ) 2 + 2J 2 (A\ + A 2 2 )+ (6) 
+2(A 2 1 -At)(sl-e 2 2 )-2(e 2 1 e 2 2 + s 2 1 J 2 +e 2 2 J 2 )] = 

can be solved in radicals, e.g. by Ferrari's method. The 
most interesting with relation to the quantum mechanics 
is the symmetric case of unbiased qubits (ei =£2 = 0) with 
maximal superposition and entanglement effects. In this 
case Eq. ^ is simplified to a biquadratic equation and 
the two-qubit system has the eigenenergies 



£ s : = tV^ 2 + (Ai±A 2 ) 2 , i = 1,2,3,4. (7) 

We focus on the still more symmetric and indicative case 
of unbiased qubits with equal tunnel amplitudes (Ai = 
A 2 = A = AZ?oi/2), in which the eigenvalues of H\ 2 arc 



E^-yJj 2 +4A 2 , E 2 = -J, E 3 = J, E^yJ.1 2 +4A 2 , 

(8) 

and the corresponding normed eigenstates are 

|*i)=o|tt) + 6|t4.> + &l4.t)+o|44>, 2(a 2 + 6 2 ) = l, 

|* a > = -±= (|n) - in)) ; 1*3) - ^= (Itt) - III)) ; 

\* i )=a'\n)+b'\n)+b'm) + a'\U), 2{a' 2 + b' 2 ) = 1, 

1 1 

a ' = VJ'^A* ( 



- ' h>- 1 1 , J 

\7,/ i +4A i I ' 2 l x x /J^+4A i 



The density operator of the system in pure (at zero tem- 
perature) eigenstates |^j) (i=l,2,3,4) have the form 



P=l*iX*i| 



(10) 
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For definiteness sake, we focus on the ground state 
of the system. Taking into account Eqs. ©, we obtain 
the density matrix pki = (k\p\l) (in the basis (j4}) in the 
ground state 



Pki 



/a 2 ab ab a 2 

ab b 2 b 2 ab 

ab b 2 b 2 ab 

\a 2 ab ab a 2 



(11) 



(note that components of pki are real as the initial tun- 
neling amplitudes Ai, A2 were chosen real, without loss 
of generality for considering stationary states and entan- 
glement measures). The density matrix of a pure state, 
such as (fTTj). has the only nonzero eigenvalue A = 1 (as 
p 2 = p for it and thus A 2 = A). The reduced density op- 
erators for either of the two qubits are defined as partial 
traces of p over basis vectors of the other qubit, that is 



h=Tt &) p=(k\p\k), \k) = |t) 2 ,|+) a , 
p 2 = Tr {1) p=(k\p\k), \k) = |t)i,|4)i- 



(12) 




Figure 4: The entanglement measures of the two-qubit system 
in the ground state l^i) as functions of the dimensionless 
interaction parameter Jq: l — £o(Jo), 2 — £d(Jo), 3—£s(Jo)- 



So, the reduced density matrices of the first and the sec- 
ond qubits in the system described by pip arc 

A/2 2ab\ 1 1 

Pi=P2=[ 2ab 1/2 J; X 1 = 2ab+-,X 2 = -2ab+- r 

(13) 

with Ai, A2 being the eigenvalues of the equal matrices. 

Now we proceed to analysis of entanglement measures. 
A two-qubit state 

\9)=A\tf) + B\U)+C\tf) +D\U) 

is entangled, that is by definition not decomposable as a 
tensor product of states of two qubits 

Wt) 1 + /m>i)®(7rt} 2 + <m> 2 ), 

if and only if the inequality AD — BC ^ is fulfilled, 
which follows straight from the definition (U) of the two- 
qubit basis. Entanglement measure of two-qubit states 



£0 = 4\AD-BC\ 2 , 



(14) 



is called the tangle, having such plain structure 2 ^. For 
the ground state given by ©, we have: 



£0 = 4(a^ - b 2 ) 2 = 



J 



72 



J 2 + 4A 2 1 + Jq 



2 1 



pU) 



where Jo = J/(2A) = J/AEqi is the two-qubit interac- 
tion energy J normed to the tunnel splitting AEqi of the 
qubits. 

There is a widespread density-matrix formalism mea- 
sure of entanglemet 2 ^: 



£ D = -Tv{(p-p 1 ®p 2 ) 2 }. 



(15) 



It quantifies the Hilbert-Schmidt distance between a den- 
sity matrix and a tensor-product matrix from reduced 
density matrices of the two qubits. Substituting matri- 
ces pi}, pH]) into dTSJ we get for the state |*i): 



tn = l-— 



,_> ..^4 , 256 fr8 = g o( 2 + g o) 



128, 

IT 



256, 

T 



G3) 

The relation between £ d and £q in p5ip is strict for the 
pure states 2 ^. 

One more popular measure is information-theoretic 
one called the entanglement of formation, defined as 2 ^: 

2 

£s = S (pi) = S(p 2 ) = -Tr{pilog 2 pi } = - ^2 A i lo S2 X* ■ 

i=l 

(16) 

With Ai, A 2 from (JT3J), Eq.pT>) gives £ s for the state 
The measure £s is related to the tangle £0 for the pure 
states as well 22 . As follows from p6p . the pure states 
themselves have £s = 0, seeing that their density matrices 
have one nonzero eigenvalue A = 1. On the other hand, 
the reduced density matrices of two subsystems of a sys- 
tem in an entangle pure state correspond to the mixed 
states. This property of efficient states of the subsys- 
tems is treated as their mutual "measuring" each other 
as parts of an entangled system. 

All the entanglement measures £q, £d and £s vanish in 
the tensor-product states, being positive otherwise, and 
maximize to unity in the most entangled states (such as 
states of the Bell's basis). In Fig.0] the £o,£d,£s a s 
functions of the qubits' interaction parameter Jo, given 
by Eqs. (JUJ, Qjft, pH), are shown for the considered 
two-qubit system in the ground state Due to sym- 

metry, exactly the same expressions hold for the highest 
excited state |* 4 ). At J o = (o=6=l/2, a' = -6' = l/2), 
i.e. at vanishing interaction strength between the qubits, 
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the states |\E'i) 
form 



| $4) decompose into the product of the 



l*i> = ^(|t)i + |4)i)®^(|t) a + |4) a ), 

I*4) = ^(|t>l-Wl)®^(lt) 2 -Il> 2 ), 



(17) 



so that all the entanglement measures vanish, the qubits 
are disentangled. Note, the reduced density matrices p\ 
and p 2 in Eq. (fT3)l describe pure states in this case (oth- 
erwise mixed ones). With increasing absolute value of 
Jo, all the entanglement measures steadily increase. In 
the characteristic point J = 1: £o(l) =0.5, <?d(1) =0.4, 
£s(l)=0.6. At J ->oo, £ = l-Jo 2 +o(Jo A ). And so, 
at increasing | J|^>A (FM-coupling), the ground and the 
highest excited states tend to 



l*i> 



-^(ltt)+iu)), i* 4 >^^(in>+nt», (is) 



the Bell's states having the maximal unity entanglement 
(while the eigenstates 1^2), 1^3), being the Bell's states, 
have the unity entanglement at arbitrary J^O). 

The eigenstates 1^,) §§§ are stationary, varying in 
time as exp(—iEit/h) with the eigenenergies Ei ^ 
(i=l,2,3,4). So the entanglement of the eigenstates is 
constant in time, as seen from Eq.(fT4"|) for Let's con- 
sider a simple dynamical model in which the two-qubit in- 
teraction under study can be turned on and off instantly. 
It simulates a situation, when interaction turning on/off 
timespan is much less than the system's characteristic 
dynamical time. Suppose, the interaction J is turned 
on between the qubits which are initially disentangled, 
specifically, in the product state from Eq. (fT7"|) . Ap- 
plying eigendecomposition of the Hamiltonian, we get the 
time- varying state vector 

mt))=e- iA ^ h \m) , (k W)) = (4 B, B, A), 

A = a(a + b)e- iE ^ h + a\a' + b')e lE ^' h , (19) 
B = b(a + b)e-' lElt/h + b'(a' + b')e lElt/h , 



and the tangle 



J2 



I 2 



'(cos(27rf/^o) - l) 2 



1 



I 2 



sin 2 (2nt j to) 



^ (20) 

where t — h/(2y/ J 2 + AE^ ) is its oscillation period. 
At Jo < 1, time of transition from the initial dis- 
entangled state to the maximally entangled one, with 
£>™« = [2/(J + Jq 1 )} 2 , is t E = t /2 (otherwise, t e = 
(t /27r)arccos(-J ( 7 2 ) > t /4 and £^ ax = 1). So, from 
a dynamical viewpoint, characteristic time of entangling 
of two initially disentangled subsystems is determined by 
their interaction energy J: t e ~ H/J at J ~ A. And 
the mere fact of entanglement between subsystems be- 
comes possible because of their interaction (£o(t) = 
at J = 0). Suppose that the interaction Jo = 1 is 
turned off at the moment t e . Then the state vector 
(k |\&(t e )) = (0, -A=, 7^)0) will remain constant in time, 
with £{){t) = 1, that illustrates the conception of quantum 



nonlocality (though the latter suggesting that consis- 
tent non-paradoxical description of interaction of quan- 
tum systems requires the theoretical setting of relativistic 
quantum field theory^). 



IV. SUMMARY 

We have analyzed the quantum inductance coupler 
based on a superconducting loop with ScS quantum point 
contact, which is intended to provide the tunable ZZ in- 
teraction between superconducting flux qubits. The fol- 
lowing features of this coupler should be pointed out that 
favorably distinguish it from the analogues by principle 
of operation. These are (i) the relatively small opera- 
tion range A/ e ~ 0.01 of the coupler controlling flux, 
ensuring the wide range of strengths J(/ e ) and hence 
promoting reduction in operational times of the coupler; 

(ii) almost symmetric form of the function J(/ e ) rela- 
tive to the point of its vanishing in the operation range; 

(iii) large attainable absolute values of the intcr-qubit 
interaction strength | J(/ e )/fcs| ~ IK. These features of 
the QUINC with ScS QPC make it optimal for coupling 
the superconducting niobium flux qubits with ScS QPCs 
as well having the tunnel splittings AEoi/^b ~ 1 — 2K, 
or AEoi/ft ~ 20-40GHzi£ (notice that properly de- 
signed superconducting Josephson circuits with QPCs 
may serve as high-quality qubits, qutrits, couplers and 
quantum detectors^ thus having great potential for the 
quantum information science). In such a way, having 
available flux qubits with the splittings AEoi/fce — 
and couplers providing energies |J//cb| — IK, a multi- 
qubit quantum system with tunable inter-qubit coupling 
energies can be constructed, which will behave strongly 
coherently at temperatures T~ 10 _2 K, enabling imple- 
mentation of quantum gates. 

We showed behavior of the widely-used entanglement 
measures for the eigenstates of the symmetric two-qubit 
system with ZZ inter-qubit interaction as functions of the 
interaction strength J. With different formal definitions, 
all the considered measures behave qualitatively simi- 
larly in agreement with physical intuition, viz. steadily 
increase with increasing J in the system's ground and 
highest excited states. Also time variation of the entan- 
glement (^-measure) in the analyzed two-qubit system 
was illustrated within an elementary model of piecewise- 
constant J(t) function. This instance shows that time- 
controlling of the entanglement is not simple problem re- 
quiring modeling of the time- varying dependence of inter- 
qubit interaction strength and realizing it experimentally 
to perform desired two-qubit gates. 
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